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Abstract. An extension of the Doi-Peliti formalism for stochastic chemical kinetics 
is proposed. Using the extension, path-integral expressions consistent with previous 
studies are obtained. In addition, the extended formalism is naturally connected to 
orthogonal polynomials. We show that two different orthogonal polynomials, i.e., 
Charlier polynomials and Hermite polynomials, can be used to express the Doi-Pcliti 
formalism explicitly 
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1. Introduction 

Stochastic chemical kinetics has been widely used in various research areas. For example, 
react ion- diffusion systems based on the stochastic chemical kinetics have been used in 
nonequilibrium physics pQ. In addition, recent development of experimental techniques 
enables us to observe various reactions in a cell in detail. It has been clarified that 
some chemical reactions in cells should be treated as systems with discrete states, and 
continuous approximation cannot be used [2]. Hence, from the theoretical viewpoint, it 
is important to develop tractable analytical schemes to deal with the stochastic chemical 
kinetics. The analytical schemes will give deep understanding for the nonequilibrium 
phenomena, and it is expected that such analytical schemes lead to rapid computational 
method for parameter estimations, which is useful for the analysis of experimental data. 

The Doi-Peliti formalism is known as a useful method to treat the stochastic 
chemical kinetics [3H5]. In the formalism, bosonic creation and annihilation operators 
are used, and we can employ an analogy between the Doi-Peliti formalism and quantum 
mechanics. Actually, there are many researches for the Doi-Peliti formalism, in which 
various techniques developed in quantum mechanics are used in order to investigate 
classical stochastic process; we can use perturbation calculations [6], renormalization 
group analysis [lj, and system-size expansion [7]. 

Although the usefulness of the Doi-Peliti formalism has been shown in many studies, 
there are some unclear points of the formalism. For example, is it possible to write state 
vectors, which are used in the Doi-Peliti formalism, in terms of analytical functions? 
If so, how should we interpret creation and annihilation operators? While it has been 
pointed out that the Doi-Peliti formalism is equivalent to generating function approach 
or Poisson representation [8], the correspondence would not be unique. If we have various 
concrete expressions for the Doi-Peliti formalism, it is possible to choose an adequate 
one depending on one's objectives; an expression may be suitable for perturbation 
calculations, and another may be tractable for numerical calculations. 

In the present paper, we firstly propose an extension of the Doi-Peliti formalism. 
While the extension includes an additional parameter, the introduction of the parameter 
does not change the coherent-state path-integral formula; all analytical techniques in 
previous works are available. The extension affect only the concrete expressions for 
state vectors in the formalism. Secondly, we point out that the one-parameter extension 
of the Doi-Peliti formalism is naturally connected to orthogonal polynomials. We find 
that two different orthogonal polynomials, the Hermite polynomials and the Charlier 
polynomials, are used to describe the state vectors in the formalism. 

This paper is organized as follows. In section [2} a brief review of the Doi-Peliti 
formalism is given. The correspondence between the Doi-Peliti formalism and the 
generating function approach is also explained. In section [3l we give the first main 
result in the present paper, i.e., the one-parameter extension. Section H] is the second 
main result; the Hermite polynomials and the Charlier polynomials are explained, and 
concrete expressions of creation and annihilation operators are given. Section [5] gives 
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discussions and concluding remarks. 

2. A brief review of previous studies 

2.1. The Doi-Peliti formalism 

In the Doi-Peliti formalism, the following bosonic creation operator a) and annihilation 
operators a are used: 

[a, a'] = aa) — a) a = 1, [a, a] = [a , a'] = 0, , (1) 

where [-, •] is the commutator, and the actions of the creation and annihilation operators 
for ket vectors \n) are defined as 

a)\n) — \n + 1), a\n) = n\n — 1). (2) 

Here, the vacuum state |0) is characterized by a\0) = 0. ([1]) means that the creation 
operator and annihilation operator do not commute. While the actions of two operators 
on the ket vectors are defined as (121), actions on bra vectors (n\ are defined as follows: 



(n\a = (n + 1|, (n\a* — {n — l\n. (3) 

The inner product for the bra and ket vectors is given as 

(m\n) = n\ 6 m>n . (4) 

where 5 m , n is the Kronecker delta. 

In order to explain the usefulness of the Doi-Peliti formalism, we use a simple birth- 
coagulation process as an example. The birth reaction X — > X + X occurs with the 
rate constant a for each particle, and its backward reaction (coagulation) X + X — > X 
occurs with 0. The master equation for the process is written as 
d 

— P(n, t) = a[(n - l)P(n - 1, t) - nP(n, t)] 

+ /3[n(n + l)P(n + 1, t) - n{n - l)P(n, t)}, (5) 

where P(n, t) is a probability of finding n particles at time t, and we define P(—l,t) = 0. 
The remarkable idea of the Doi-Peliti formalism is the usage of a single vector \ip(t)) 
which is a collection of a series of infinite number of P(n,t): 

CO 

W*)> = j>(n,t)|n), (6) 

n=0 

Using the vector \ip(t)), the master equation fl5]) is rewritten in a compact form 

^ t \m) = L(a\a)\m), (7) 

where 

L(a\ a) = ot{a) - l)o)a + (3(1 - a))a ] a 2 . (8) 

Because of the similarity with quantum mechanics, the Doi-Peliti formalism is also 
called the second-quantization method or the field-theoretic approach. In addition, 
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similarities between the Doi-Peliti formalism and quantum mechanics enable us to use 
various analytical techniques in quantum mechanics in order to investigate classical 
stochastic processes, as explained in section CD 

Here, we briefly explain some basic definitions for a coherent-state path-integral 
expression, which is especially useful in the Doi-Peliti formalism. In order to derive the 
path-integral expression, coherent states and a decomposition of unity play essential 
roles. The coherent states are defined as 



oo 1 

n=0 
oo 1 



n\[z 



(9) 



(10) 



71=0 



where z is a complex number, and z* is the complex conjugate of z. Using the coherent 
states, the decomposition of unity is obtained as follows: 

^ ' d 2 z 



1 = y^-r|ra)H = — -.\n){m\8 n , m = \ 
z — ' nl A — ' A — 4 ml I 

n=0 n=0 m=0 J 



7T 



\z)(z\ , 



where we used 



d 2 z 



-\z\ 2 z * m z n 



(11) 



(12) 



with the integration measure d 2 z = d(Re z)d(Im z). Using the decomposition of unity, 
it is straightforward to obtain the path- integral expression; for details, see Ref. [I]. 

Note that there are some differences between the Doi-Peliti method and quantum 
mechanics. One of them is the calculation scheme for expectation values of observables. 
Different from quantum mechanics, an expectation value of observables in the Doi-Peliti 
formalism is obtained by using a projection state: 



oo 1 

<^<ok = ]T-H. (is) 

n=0 

For example, the average of n is given by Yl^Lo n>P(n, t) = (V\a^a\ip(t)). 
2.2. Correspondence with the generating function approach 

While the Doi-Peliti formalism is one of the analytical methods to treat stochastic 
chemical kinetics, the generating function approach is a more well-known method and 
has been studies a lot 0. For the previous example in section I2.1[ the generating 
function G(x, t) is defined as 



G(x,t) = J2 p (n,t) 



x 



n=0 



and the time evolution equation for G(x, t) is written as 



dt 



G(x,t) 



a(x — 1)^-7- + 0(1 — x)x-^r 
dx dx 2 



G{x,t). 



(14) 



(15) 
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The above partial differential equation can be used instead of the simultaneous 
differential equations (jSJ). 

One may consider that ( fT5l) is similar to (J7|) with (JHJ). In fact, as written in 
section [Q, it has been pointed out that there is a correspondence between the Doi-Peliti 
formalism and the generating function approach [8J. When we interpret the creation 
and annihilation operators as 

a) = x, a = — , (16) 
ax 

the linear operator L(a),a) for the time-evolution of \ip(t)) in the Doi-Peliti formalism 
is obtained from that of G(x, t) by replacing x and 4- with and a, respectively. In 
addition, explicit expressions for the ket and bra vectors in the Doi-Peliti formalism are 
obtained as follows: 

dx5{x)( — j (•)• (17) 
This correspondence has been used to discuss duality relations in stochastic processes 

nm. 



3. One-parameter extension of the Doi-Peliti formalism 

Starting from the same commutation relation ([I]) and the same definitions for the ket 
vectors ([2]), it is possible to add one parameter A to the formalism; this extension is one 
of the main results in the present paper. That is, we define the following inner product 
instead of @: 

(m\n) = X n n\ 5 m>n . (18) 

According to the replacement of (Hj) with ([TBI , the actions of the creation and 
annihilation operators on bra vectors (n| change as follows: 

(n\a = in + 1|A~\ (n\o) = (n - l|nA. (19) 

The derivation of ( Tl9|) is as follows: Firstly, we have {n + l\a^\n) = (n + l\n + 1) = 
{n + l)A n+1 because o)\n) = \n + 1). Secondly, we assume the action of the creation 
operator on the bra vector as (n + l|a' = (n\a, where a is a scalar value. Then, 
(n + l|a^|T2,) = a(n\n) = an\X n , and we have a = (n + 1)A. Hence, the second equality 
in f JT9|) is obtained. Using the similar discussions, the first equality in ([19]) is easily 
checked. 

We here note that the path-integral expression must not be changed due to the 
above one-parameter extension because the final expression of the path-integrals consists 
of integrals only for parameters in the coherent states, i.e., z and z*\ the final expression 
does not depend on the definition of \n) and (n\. Actually, there is no need to change 
the definitions of the projection state and the coherent states. Since the actions of the 
creation and annihilation operators for the bra vectors are modified as ( Tl9|) . we have 
\n) = (a^) n \0) and (n\ = (0|(aA) ra . Hence, from the same definitions with the usual 
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Doi-Peliti method, we obtain slightly different expressions for the projection state and 
the coherent states when we write them explicitly using the bra vectors (n\ as follows: 

oo 

n=0 

oo 

|z>=e" t |0> = V-z»|n>, (21) 
' m 

n=0 

W^^-E^H U) ■ (22) 

n=0 ' v 7 

The above discussions suggest that there is no need to change the definitions of the 
projection states and coherent states. In addition, the decomposition of unity is 
calculated as 

00 1 1 r a% 

l = £^>M=/^e-'*'V)<4 (23) 

71=0 ' J 

The unity in the extended Doi-Peliti method has the same expression with the usual 
Doi-Peliti method in terms of the coherent states, and therefore we obtain the same 
path-integral expressions even in the one-parameter extension, as expected. 



4. Correspondence with orthogonal polynomials 

As explained in section 12.21 the operators and vectors in the Doi-Peliti formalism are 
explicitly rewritten in terms of polynomials and differential operators. The expression is 
not unique, and we here show two different expressions for the one-parameter extension 
of the Doi-Peliti formalism. 

In the following discussions, we restrict the additional parameter A as a positive 
real variable, i.e., A > 0, in order to see the connection to the orthogonal polynomials. 



4-1. Hermite polynomials 

One of the expressions is obtained from the Hermite polynomials [llj. The Hermite 
polynomials are defined as 

Jn 

H n (x) = {-l) n e* — e-* 2 , (24) 

where n G N, x G K. Introducing a scaling variable A G R, we define the following 
rescaled Hermite polynomials: 

m*) - /gf B. (^) . (25) 

Using the property of the Hermite polynomials, it is straightforward to verify the 
following three-term recurrence formula: 

H^(x) = xH£\x) - XnH^x). (26) 
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In addition, the rescaled Hermite polynomials satisfy the following orthogonality 
relation: 

/+oo 
m n x \x)H^(x)^ x \x)dx = X n n\5 ntm , (27) 
-oo 

where 

/i (A) (x) = _L e - 2 /(2A)_ (2g) 

V^x 

As one can easily see, the orthogonality relation f[2"T|) corresponds to the inner 
product (I18p in the one-parameter extension of the Doi-Peliti method. Actually, if we 
define 

/oo 
dx^ x \x)Hi x \x), (29) 
-oo 

a) = x — A — , a = -J— • (30) 
ax ax 

all properties in the one-parameter extension of the Doi-Peliti formalism are recovered. 
For example, the action of the creation operators on the bra vector, (n\a* = (n — l|nA, 
is verified by using the recurrence formula ( 1261) and a partial integral. 



4-2. Charlier polynomials 

Another expression is obtained from the Charlier polynomials [llj. The definition of 
the monic Charlier polynomials is 



n , 

ci A) (^)=E(- A ) n "* a?( * ) ( 

k=0 ^ 



(31) 



where x^ = x(x — 1) • • ■ (x — k + 1) and n E N, x G N. Note that the variable x is not a 
real value but a natural number, which is different from the Hermite polynomials. The 
Charlier polynomials satisfy the recurrence formula 

Cl%(x) = (x-n- X)C^(x) - XnC^x) (32) 

and the orthogonality relation 

E C4 A >(s)CW(s)^e-* = \n n \5 m>n . (33) 
L — ' x\ 

x=0 

We here introduce the following definitions for the bra and ket vectors: 

00 \ x 

\n) = C^(x), H^-e- A Cf)(:r), (34) 

In order to recover the properties of the one-parameter extension of the Doi-Peliti 
formalism, we define the creation and annihilation operators as 

a^f(x) = xf{x-l)-Xf{x), af(x) = f{x + l)-f(x). (35) 

Some techniques in discrete mathematics is needed to verify (JT9l) for the above 
definitions, and the explanations are a little complicated. We will give them in 
Appendix. 
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5. Discussions and concluding remarks 

We have found a one-parameter extension of the Doi-Peliti formalism, and the extended 
formalism is deeply related to orthogonal polynomials. Although the correspondence 
with the generating function approach has already been known previously, essentially 
different expressions for the one-parameter extension of the Doi-Peliti formalism have 
been obtained. 

While the extension does not affect the coherent path-integral formulation, the 
additional parameter will be useful for some cases, for the following reason. The time- 
evolution operator L(a^,a) for the state vector \i[)(t)) is determined from the master 
equation, and it is independent from the additional parameter A. On the other hand, 
the parameter A can change analytical expressions for the bra and ket vectors, and hence 
analytical expressions of the state vector \ifj(t)) depend on the parameter A. As a result, 
it is expected that one can choose an adequate parameter A which is suitable to express 
the solution \ijj(t)). Actually, it has been known that some stochastic processes are 
related to orthogonal polynomials, and the choice of the parameters in the orthogonal 
polynomials is important to express the analytical solutions [12J. 

In the correspondence between the generating function approach and the usual Doi- 
Peliti method, the expression with the Dirac's delta function has been known. However, 
the delta function may be intractable when we want to construct some numerical 
methods based on the formalism. In contrast, the Hermite and Charlier polynomials 
would be tractable, and hence they will become helpful in numerical computations. As 
explained in section [H the Doi-Peliti formalism has been widely used. We expect that 
the correspondence to the orthogonal polynomials will give useful methods in order to 
investigate nonequilibrium behaviour. 
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Appendix A. Validity of the definitions in section 14.21 

It is well-known that the difference operator [13] 



so that a\n) = n\n — 1) is verified. In addition, combination of the recurrence formula 
(152"]) and the difference equation [TT] 



Au(x) = u(x + 1) — u{x) 
acts on the Charlier polynomials as 



(A.l) 



AC?W(x) = Cj»(x + 1) - Cf )(*) = nC^ix), 



(A.2) 



— n 



CW(x) = AC< A) (i + l)-(i + fl 



)Ci x \x)+x 



(A.3) 
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gives 



ci X } 1 (x)=xCi x \x-l)-XC^(x), 



(A.4) 



which corresponds to a*\n) = \n + 1). 

In order to check (n|o' = (n — l\n\, a partial summation is available. Using the 
shift operator defined as 



Eu(x) = u(x + 1), 
the partial summation is given by [13] 

uAu = uv — A-uEf. 
We here note that we can interpret a'\n) as 

^(rc - 1) - XC { n x \x) = -xAC {x \x - 1) + xC^^) - XC { n x \x). 
Hence, (n\o)\m) is expressed as 

oo _^ . x 

E ^T^^) h^ffC* - 1) + *Cl A) (z) - ACff(: 

z=0 



(A.5) 
(A.6) 



x! 

oo 



-C<»(x)(-x)C<»(x-l\ 



x=0 



x=0 L 



A 



X\x 



e A 



.r! 



e- x X x 



-Ci x \x)x \ ECS?{x - 1) + -^rCi x \x) {xCW(x) - ACW(x)} 
e" A A a; 



x=0 

e~ x X x 

+ r- 

x\ 



XI 



x- 1)! 

c£\x){xC£\x)-XC£\x)} 



C {X) {x)\c£\x 



Y. e ^ X {C { n X) (x+l)-CL X \x)}C( X \x) 



x=0 
oo _^ , x 

n X! 



(A.7) 



where we used flA.2j) to obtain the final equality. Hence, the action of the creation 
operator a) on the bra vector, (n\a) = (n — l\nX, is verified. The action of the 
annihilation operator a on the bra vector, (n\a = (n + 1|A _1 can be checked in a similar 
manner. 



References 

[1] Timber U C, Howard M and Vollmayr-Lee B P 2005 J. Phys. A: Math.Gen. 38 R79 

[2] Rao C V, Wolf D M and Arkin A P 2002 Nature 420 231 

[3] Doi M 1976 J. Phys. A: Math. Gen. 9 1465 

[4] Doi M 1976 J. Phys. A: Math. Gen. 9 1479 

[5] Peliti L 1985 J. Physique 46 1469 

[6] Dickman R and Vidigal R 2003 Braz. J. Phys. 33 73 



One-parameter extension of the Doi-Peliti formalism and relation with orthogonal polynomialslO 

[7] Itakura K, Ohkubo J and Sasa S-i 2010 J. Phys. A: Math. Theor. 43 125001 

[8] Droz M and McKane A 1994 J. Phys. A: Math.Gen. 27 L467 

[9] Gardiner C W 2004 Handbook of Stochastic Methods 3rd edn (Berlin: Springer) 

[10] Ohkubo J 2010 J. Stat. Phys. 139 454 

[11] Chihara T S 1978 An introduction to orthogonal polynomials (New York: Gordon and Breach) 

[12] Schoutens W 2000 Stochastic processes and orthogonal polynomials (New York: Springer- Verlag) 

[13] Graham R, Knuth D E and Patashnik O 1994 Concrete mathematics: A foundation for computer 
science (Reading, MA: Addison- Wesley) 



